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Abstract 



k> , We consider the infrared (IR) limit of the nonhnear integral equations (NLIEs) for 

the boundary supersymmetric sine-Gordon (BSSG) model, previously obtained from 
the NLIEs for the inhomogeneous open spin-1 XXZ quantum spin chain with general 
integrable boundary terms, for values of the boundary parameters which satisfy a 
certain constraint. In particular, we compute the boundary S matrix and determine 
the "lattice - IR" relation for the BSSG parameters. 
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1 Introduction 

One-dimensional spin systems such as the open XXZ quantum spin chain, e.g., see [T|-[T9] 
and references herein, have been crucial in the study of ground and excited states of the cor- 
responding field theoretical models i.e., the boundary sine-Gordon (BSG) [20] and boundary 
supersymmetric sine-Gordon (BSSG) models [2lll22]- In particular, studies involving the cru- 
cial ultraviolet (UV) limit and infrared (IR) limit which include the computation of boundary 
S matrices (reflection factors) for soliton and breathers, have attracted much interests over 
the years. One effective tool for such studies is the method of nonlinear integral equations 
(NLIEs). The NLIEs have been obtained for SG model from the Bethe ansatz equations of 
spin-1/2 XXZ spin chain. These NLIEs have revealed important results about these mod- 
els. The studies initiated for periodic boundary conditions [23]- [25] were also extended to 
Dirichlet and more general boundary conditions [26] -[IQ]. They have been used to compute 
quantities such as the 5* matrices (both bulk and boundary), bulk and boundary energies, 
central charges and conformal dimensions. 

The supersymmetric counterpart, the SSG model [ll]-[5l], have also received equally 
extensive attention. Initially, a set of NLIEs for the close spin-1 XXZ chain, which is related 
to the SSG model, were derived in [551 ES] from the T-Q equations of close spin-1 XXZ spin 
chain, by exploiting the analyticity of the transfer matrix eigenvalues described by these 
equations. Such a method have been later used to derive the desired NLIEs describing SSG 
models with periodic boundary conditions [52], the Dirichlet [53] and also for a more general 
boundary condition [51] . These NLIEs were used to calculate finite-size properties of ground 
and excited states. 

Our motivation for the present work comes from these efforts. Utilizing solutions for the 
open spin-1 XXZ chain with general integrable boundary terms (rather than the diagonal 
boundary terms used in [53]) [H] that have been recently proposed [T5fH . a corresponding set 
of NLIEs was derived for a BSSG model [5l]. We thus plan to analyze the IR limit of these 
NLIEs along the line of [53] and obtain the corresponding boundary 5" matrix. We will show 
that the result agrees with that obtained by Bajnok et al. in [50]. We shall consider only the 
one-hole state in this paper. Another important objective of this work is to determine the 
relation between the boundary "lattice-IR" parameters for the BSSG case g (for other than 
that of the Dirichlet boundary condition case). Such a relation is well-known for the BSG 
case, and has been found for BSSG models with Dirichlet boundary conditions [53], labelled 
there as the Dirichlet BSSG"*" model. 



^The lattice boundary parameters satisfy a pair of real constraint. 

^Note that in |50], two BSSG models are discussed: BSSG^ and BSSG+. Our results arc related to the 
latter model by Bajnok et al. 



The outline of the article is as follows. In section 2, we briefly review the open spin-1 
XXZ quantum spin chain, namely the Hamiltonian, the T-Q equations and the NLIEs of 
the model. The BSSG model and the set of NLIEs that describe the model are reviewed in 
section 3, which are reproduced from |5l]. In section 4, we give the main results of the paper. 
We consider the IR limit of the NLIEs for a state with one hole which is subsequently used 
to compute the boundary S matrix. We then determine the relation between the boundary 
"lattice-IR" parameters. Finally, a brief discussion of our results and some open problems 
conclude the paper in section 5. 

2 Open spin-1 XXZ chain 

Since the set of NLIEs that describe BSSG models and the NLIEs for the inhomogeneous open 
spin-1 XXZ quantum spin chains are closely related, in this section, we briefly review some 
crucial results on the open spin-1 XXZ quantum spin chain which include the Hamiltonian 
of the open spin-1 XXZ chain, the T-Q equations for the inhomogeneous model, important 
auxiliary functions and the NLIEs. The readers are urged to refer to ^5\ El] for more details 
on these materials. 

2.1 Hamiltonian and T-Q equations 

We begin with the Hamiltonian of the open spin-1 XXZ chain. We adopt the notations used 
in [15], 

N-l 

n = Y,Hn,n+l + Hb, (2.1) 

n=l 

where Hn^n+i represents the bulk terms. Explicitly, these terms are given by [57] . 
Hn,n+i = cr„ - M' + 2 sh^ [< + {S:,f + {S:^,f - «)2] 

- 4sh2(|)K< + aX^) , (2.2) 



where 



O'n — bji ■ ^n+1 1 0"„ — ^n^n+1 + ^n^n+1 ) ^n ~ ^n^n+1 i (2-3) 



and S are the su{2) spin-1 generators. Hb represents the boundary terms which have the 
following form (see e.g.. 



+ arS'^ St + agS^ S'^ + {a^ ^ bj and I ^ N) , (2.4) 



where S^ = S^ ± iS^. The coefficients {oj} of the boundary terms at site 1 are functions of 
the boundary parameters (a_, /3_, 6'_) and the bulk anisotropy parameter rj. They are given 
by, 



as 

a4 

as 



-ao (ch 2a_ — ch 2/3_ + ch r^) sh 27] sh 77 . 



-oq sh 2a_ sh 2/3_ sh 27] , 



.2e_ 



— -aoe "sh2?7sh?7, 
8 



1 



— aoe " sh2r7sh?7 
8 

aoe^~ 

—aoe' 
—aoe 



77 77"\ 3 

ch /3- sh a_ ch — h ch a_ sh /3_ sh — j sh 77 ch 2 77 , 

?7 77 \ 3 

ch /3_ sh a_ ch — h ch a_ sh /3_ sh — j sh 77 ch2 t] , 

77 77 \ 3 

ch /3_ sh a_ ch ch q;_ sh /3_ sh - j sh 7/ cha 7/ , 



77 ?7 ' - 

ch /3_ sh a_ ch ch a_ sh /3_ sh — ) sh 77 ch 2 7/ . 



where 



ao 



^^ 



r/. 



V^ 



T]. 



sh(a_ - ^) sh(a_ + j) ch(/3_ - ^) ch(/3_ + ^) 



(2.5) 



(2.6) 



Similarly, the coefficients {6j} of the boundary terms at site A^ which are functions of the 
boundary parameters (a+,/3+,6'+) and r/, are given by the following correspondence, 



a_ — >a+ ,/3_ ->— /J^ 



>e+ 



(2.7) 



We stress beforehand that the NLIEs in [51] were derived from the T-Q equations for 
the (inhomogeneous) open spin-1 XXZ chain with integrable boundary terms as given in 
[T5] . namely when the boundary parameters (ct± , /3± , 6'±) obey the following constraint 

I9i[in],[i5]i 



a_ + /3_ + a+ + /3+ = ±{d^ - 9+) + rj . 



{21 



As adopted in [38] , a convenient redefinition of bulk and boundary parameters as follows 



7] = Ifi . 



a± = iiia± 



/3± = nb± 



e 



± = l^C± 



•^More generally, the last term in (|2.8p is krj, where k is an odd integer. Following 
See the reference for details. Also refer to [15] for more complete discussion on this. 



(2.9) 
], we take k = 1. 



results in the following pair of real constraints: 

a- + a+ = ±|c_ — c+l + 1, 

b_+b+ = 0. (2.10) 

where fJ. ,a± ,b± ,c± are all real, with < ;U < |. 

Next, we review the model's T-Q equations that describe the eigenvalues of the commut- 
ing transfer matrices [15]. Two relevant commuting transfer matrices for spin-1 XXZ chain 
are Ti[u) with a spin-i (two-dimensional) auxiliary space, and T2{u) with a spin-1 (three- 
dimensional) auxiliary space. The eigenvalues of Ti{u), Ai(m), can be written as (following 
[53] and adopting the notations used there) 

Ai(m) = /i(m) + /2(m), 

h(u) = sinh(2M + t^)B^+\u)(l)(u + t^) ^^"!''^^ , 

Q{u) 

hiu) = sinh(2M - i^)B^-\uUiu - 1^) ^^"^^ l^'' , (2-11) 

Q{u) 

where 

(j){u) = sinh^(M-A)sinh^(M + A), 

B^^) (u) = sinh(u ± ^J^—^) smh{u ± ^J^-^) cosh(M t ^^^-^) cosh(M t ^^^— ) , 

M 

Q{u) = 1 I sinh(u — -Ufc) sinh(M -|- Cifc) . (2.12) 

fc=i 

As in [S3] , the lattice boundary parameters have been redefined as A± = 2a±. — 1 , B±. = 
2ib± + 1. A is the inhomogeneity parameter, A^ is the number of spins and M = N represents 
the number of Bethe roots which are also the zeros Vk of Q{u). For this particular model of 
the XXZ chain, one generally needs to consider two groups of transfer matrix eigenvalues, 
labelled as A*^ 2 '!)(=•=) (y) in [TJ], to obtain all 3^ energy eigenvalues. However, we restrict to 
only one of them that contains the ground state, namely M2'^)^'(^u). Refer to [I5] for greater 
details on this. From the fusion relation, the eigenvalues of T2{u), A2{u) can be written as 
(see e.g. [T5] ) 

A^iu) = A,{u- '-^) A,{u + |) - f{u) (2.13) 

where 
f{u) = (j){u+^)(f){u-^)smh{2u-2tfi)smh{2u + 2tfi)B^+\u+'-^)B^-\u-'-^) 

= /i(n + |)Mu-|) (2.14) 



Using fl2.14p . fl2.13p can thus be written as 

= sinh(2M)A2(M) (2.15) 

where 
A2H = sinh(2n-2^/i)5(-)(n-^)5(-)(n + ^)0(M-^)0(M-^^^^^^^tM 



2' ' 2'^' 2 '^' 2' Q{u 



^/i^ 



+ sinh(2^)BW(n - '-^)B^-\u + ^)0(« - ^)0(« + !/^) ^(^ + f^) ^(^ ' f^) 

+ sinh(2« + 2^^^)m^\u + |)5W(« - |)0(« + ^)0(« + f )^^ 

= Ai(m) + A2(m) + A3(m) (2.16) 

From (I2.15P and (I2.16p . one has the following 

(2.17) 

The next crucial review are the following auxiliary functions b{u) and b{u) that enter the 
NLIEs: 

Ku) = -— , b{u) = b{-u) = ^— (2.18) 

M{u) \i{u) 

For real u, b{u) is the complex conjugate of b{u). Using (12. lip and ( I2.17p . b{u) becomes 

^.._ Ai(^-f) ^in-^^) ^(-)(n + f) g(^ + f) 

^' smh{2u + 2ifx)(P{u+'f)(P{u + ^)B(+){u-f)B(+){u+f)Q{u-^) ^' ' 

Finally, we conclude this section with another important auxiliary function, y{u) that enters 
the NLIEs, 

, , sinh(2-u)A2('u) 

fw 

2.2 NLIEs for inhomogeneous open spin-1 XXZ chain 



Next, the NLIEs are reviewed as presented in [54J. It is sufficient here to mention that the 



NLIEs are obtained utilizing the analyticity of In A2(m) = In l^V near the real axis, namely. 



0= (b du [lnA2(u)]V'=" (2.21) 



where the contour C is chosen as in figure below, e is small and positive. Refer to analysis 
in [5lj for details. 



-le 



Co 



Integration contour 



k{u) is any function whose only real root is a simple zero at the origin, that is fi;(0) = 
, «'(0) 7^ 0. The prime denotes differentiation with respect to u. These NLIEs are written 
in coordinate space: 



\n.h{u) 



du' G{u -u' - ie) ln(l + h{u' + it)) - / du' G{u - m' + ie) ln(l + b{u' - ie)) 

I sinh 



+ / du G2{u — u + ie) \n{l + y{u — ie)) + i2N tan 



+ i du' R{u') — in . 



7VU 

^^ 

[ cosh ^ 



/OO /"CXD 

du' G{u -u' - ie) ln(l + h{u' + ie)) + / du G{u -u + ie) ln(l + h{u' - ie)) 
-OO J —OO 

r°° ( sinh — 

+ / du G2[v! - u ^ ie)\Yy{\ ^ y[u ^ ie)) - ilN twy-^ \ ^ 

i-oo ycosh^^ 

r-u 

— i du' R{u') + in , 
Jo 

/OO /"OO 

du G2{u — m' + ie) ln(l + h{u' — ie)) + / du G2{u — u + ie) ln(l + b{u' + ie)) 
OO J —OO 

r-u 

+ 47ri / dv! G2(-m') , (2.22) 

where (?(«) and G2{u) are the Fourier transforms of G{k) and G2ik) respectively Q These 



Conventions used are 



iku r/^,\ Ji^, ff:,<\ I ^ — iku 



f{k) = I e^'=« fiu) du , fiu) = — / e-^'^" /(fc) dk . 



are given below, 

-, . sinh ((vr- 3u)|) 

G{k) = r^ , ^^^^ ,, , 2.23 

^ ^ 2coshf sinh((7r-2/i)|) ^ ^ 

G2{k) = J~^ ^ , (2.24) 

e 2 + e 2 

Similarly, R{u) refers to the Fourier transform of R{k) which is given below, 

s+ sinh ((vr - /i|A+|)|) + s_ sinh ((tt - /i|y4_|)|) 



R{k) = 27r< 



2cosh(f )sinh((7r-2/i)f) 



sinh (|/ig+) + sinh (|/xi?„) cosh(f)sinh((3/i-7r)|) ^ 

2cosh(f)sinh((7r-2/i)|) cosh(f ) sinh ((2/x - 7r)|) '^ ' ^' ' 



where s± = sgn(y4±) or in terms of a± is sgn(2a± — 1). 



3 The BSSG model and the NLIEs 

In this section, we briefly review the BSSG model, mainly reproduced from ^22j. In addition, 
we also give the corresponding NLIEs. 

3.1 The BSSG model 

The Euclidean-space action of the BSSG model is given by 

/oo /»0 /»oo 

dy dx Co + dy Ci,, (3.1) 

-oo J —oo J —oo 

where the bulk Lagrangian density is given by 

Co = 2dz(pdz(p — ^tpdzi/j + ^ipdzilJ — 4 cos Lp — A-ipijj cos — , (3.2) 



In fl3.2p . ip and ip are the two components of a Majorana 
z = X — iy. The boundary Lagrangian at x = is given by 



ermion field, and z = x + iy, 
Ch = ipTp + iadya — 2p{ip)a{ip — ■?/') + B{(p) , (3.3) 



^This corresponds to the BSSG"'' models in [SD], which we are interested in here. 



where a is a Hermitian Fermionic boundary degree of freedom. The functions p{<f) and 
B{(p), which are potentials that are functions of the scalar field ip, are determined from the 
requirement of boundary integrability and supersymmetry. They are found to be [22] 

1 \/F' 1 D V sin — 
B{ip) = 2v cos -{cp - ipo) , p{ip) = ^ sin -{if - D) where tan — = yo^_^ (3.4) 

where F = ^Jv'^ — 8t; cos ^ + 16 

The parameters v and y^o are arbitrary and real. In the limit that the boundary mass 
parameters tend to infinity, one arrives at the BSSG models with Dirichlet boundary condi- 
tions studied in [53] (labelled as Dirichlet BSSG+ there). Such a BSSG model corresponds 
to the open spin-1 XXZ chain with diagonal boundary terms [13J. 



3.2 NLIEs for BSSG model 

Our final review is the set of NLIEs that describes the BSSG model. In the continuum limit 



which consists of taking A — )■ oo, A^ — )• oo and lattice spacing A — t- 0, such that the interval 
length L and the soliton mass m are given by 

L = N/\, m = —e'^, (3.5) 

respectively, (as stated in earlier works [52|l53]) the NLIEs of the inhomogeneous spin-1 XXZ 
chains lead to the NLIEs of the SSG models. With the following change of variable 

e = ^. (3.6) 

(1^^ becomes 

/OO /'OO 

d9' G{e -9'- is) ln(l + hie' + ie)) - / dO' G{e -9' + le) ln(l + h{9' - %e)) 
oo J ~oo 

/oo 
de' Q2{0 -e' + ie) ln(l + y(^' - is)) + i2mL sinh 9 + i PuryiO) - ivr , 
-OO 

/OO /"OO 

de' G{e -e' - ie) ln(l + b{e' + ie)) + / de' G{e -e' + le) ln(l + h{e' - ie)) 
oo J —oo 

/oo 
de' Q2{d' -e + ie) ln(l + y(^' + ie)) - i2mL sinh e - i PbdryiO) + in , 
-OO 

/OO /"OO 

de' (k2{e -e' + ie) ln(l + h{e' - ie)) + / de' G2(^' -e + ie) ln(l + h{e' + ie)) 
oo J —oo 

+ tPyie). (3.7) 



where following definitions have been used, 



/x TT vr 



Moreover, G(^) and (j2{0) are taken to be 

TT TT 



^ , , dke-'^^"^'^ G{k), (3.9) 



oo 



G2W = ^G,{^) 

TX 71 



/^ 



^.ie-*-/-G.W = ^— ^, (3.10) 

where G{k) and G2{k) are as given in fl2.23p and fl2.24p respectively. Definitions of Fourier 
transform as given in footnote 4 have been employed. Similarly, Pbdry{(^) and Py{0) are given 
by 



bdryKP ) J^^bdry\ 

vr 



/6 /"OO 

rf^' / dke-'^^^''^ R{k), (3.11) 

-6» J-00 






and 



p.(^) = n(^; 



= Ait dO' G2{-e') = -2i\ntanh- -271, (3.12) 

J -00 2 

respectively. In (13.111) . R{k) is given by (I2.25p . As will be shown in subsequent sections, 
boundary terms Pbdryi^) and Py{9) are essential when computing the boundary S matrix. 



4 IR limit 

In this section, we consider the IR limit of the NLIEs reviewed in section 3.2, namely the 
set of equations listed in (13. 7p . The IR limit consists of taking mL — > 00. In this limit, the 
NLIEs for a one-hole state with rapidity 6h becomes, 

Inb(^) = i2mLsmhe + iPbdry{O)+iQ{O-0h) + ig{O + eh) 

/oo 
de' (k2{e -e' + le) ln(l + y{e' - is)) - m , (4.1) 

-00 

iny(^) = tPy{e) + igy{e-eh) + tgy{e + eh), (4.2) 



where g{9) and gy{0) are the hole source terms (see [52l 155] ) 

g{9) = 2n d9' G{9') , gy{e) = -ilntanh- + - . (4.3) 

Jo 2 2 

As pointed out and employed in |52l [53] for cases with periodic and Dirichlet boundary 
conditions, we stress here that when mL — )■ oo, the integral terms in (13. 7p with h{6) (and 
b{6)) are of the order 0(6""^^) and can therefore be neglected. The integrals involving y{6) 
do not enjoy such property and thus remain after the limit is taken. It is left to be shown 
in the next section that this will lead to a relation equivalent to the Yang equation for a 
particle on an interval. This will be exploited to compute the boundary S matrix. 

4.1 Boundary S Matrix and lattice-IR relation 

As for the bulk theory |171 HH], in [50], Bajnok et al. suppose that the full BSSG boundary 
S matrix is a product of the SG and RSOS boundary 5* matrices. Also refer to [53] where 
complete boundary S matrix of this form for the BSSG models with Dirichlet boundary 
conditions is given. 

In this section, we shall also find such a form here. As pointed out in [52| [55] for the 
periodic boundary conditions and in [53J for the case with Dirichlet boundary conditions, 
here we also have that In h{9h) is in times an odd integer which translates into the following 

which is obtained after (14.11) is evaluated at 6h followed by exponentiation. In (14. 4p . C{6) is 
given by 



/oo 
de' Q2{0 -0' + ie) ln(l + y(^' - is)) . 
-oo 



(4.5) 



Thus, (14. 4p should be equivalent to the Yang equation for a particle on an interval of length 
L. For an excellent discussion on this matter for the BSG model with general integrable 
boundary interactions, refer to section 4.2 of ^40j. Hence comparison of these two results 
gives the following for the product of boundary S matrices, 

R{eh ; A , r/_ , ^_) R{eh ; A , r/+ , ^+) = e^P^iry{eu)Me{2Bu)+c{eu) _ (45) 

Note that A is the IR bulk SSG parameter while ?7± and 'd± are the IR boundary SSG 
parameters. We first evaluate the factor e^^^^\ The boundary term Py{9) given by (I3.12p is 
essential in this computation. As found in [54], this coincides with the expression obtained 
in [53]. In addition, we remark that the term G2{9) (and hence its Fourier transform, G2{k)) 

10 



and equation (14. 2 p are identical to the corresponding ones in [53] for Dirichlet BSSG^ models. 
Consequently, the computation and the result for the factor e^^^^"^ are also identical and we 
therefore choose to omit the steps here. We urge the readers to refer to [521 [53] for details. 
Rather, we present only the result below, 

J^C{9f,) = ^ dk f-^ 2 , (4.7) 

d9h ^ ' 4 J_^ cosh" f cosh^ irk ^ ' 

which after integration and letting k = - yields, 

e^('''^)~P(^,)2, (4.8) 

where P{0) is given by 

'i [^ dt sin(2t^/7r) 



Indeed (14. 9p is a reflection factor of the boundary tricritical Ising model [58]- [60], the integral 
representation of which is given in [6l]. This is the RSOS factor. In [50], this factor is given 
in the form of infinite products of gamma functions (Refer to equation (14) of the reference.). 
As before in [53], this result is given only up to crossing factors of the form e^""***^. 

The remaining term e^^^'^'^y'^^^^^^^^'^^^^ in (14. 6 p will be shown below to represent the SG 
factor. At this point, we first proceed to compute -^Pbdry{Gh)- First, recalling (13. lip and 
(I2.25p . then differentiating with respect to Oh, we obtain the following 



.a-^bdry{dh) = I dke . 

^^h J-oo [ 2cosh^sinhf(^-2)f^J 



coshfsinhfc-3)f 
+ (+ ^ -) + V^ f ) ■ (4.10) 



coshf sinh('(^-2)f 



The symbol (+ f-)- — ) represents the terms with A^ — )• y4_, 5+ — )■ B_, s+ — t- s_. Similarly 
the second term gives 



I Trfc' 



J roc sinhf(^-3)^^ 

-^,{26,) = / dke-^'^'^ }' ''^ (4.11) 

d^h J-oo 2 cosh ^ sinh f (^ - 2)-^^ 



4 



11 



which together with f l4.10p results in 
d 



r .,, rksinh((7r-/.|A+|)g)+sinh(fi?+) 

Mh J-oo [ 2cosh^sinhf(^-2)^l 



2 

sinh^sinhfr^-S)^)") ^ 

+ (+«-) + 2 '- y-- ''A , (4.12) 

sinh7rA;sinhf(^-2)^j J 

Equation (I4.12p resembles the soliton reflection amplitude of the BSG model. Indeed by 
comparing with the soliton reflection amplitude P+{6 ,ri, i)) of the BSG model [20], or more 
speciflcally P+{9 , rj_^'d-)Pj^{0 , — r7+, 'd+) for two boundaries, which has the following integral 
representation, [IQ] 



1 d r°° fsinhfri-^l^l sinh f (1 + ^)^') 

i^in[p,(^,,_,^_)p,(^,-,„^,)] = / dkeM ''T^^\, :\\\K ^' ^'^'^ 



idO'"'"^'" '' '~ ''"'"' '-r^^-r/j J_^^~'~'' 1 2cosh^sinhff 2cosh^sinhff 



+ 



sinh((l + ^)f) sinh((l + ^)f 



2 cosh ^ sinh || 2 cosh ^ sinh || 



^sinh^sinh((i-l)^)) 
sinh irk sinh ■ 



+ 2 T-r-T^T^n^fc^^^ ^ (4-13) 



4A 

one will be able to determine the relation between the boundary lattice and boundary IR pa- 
rameters. In (I4.13p . A is the bulk IR parameter, ?7± and •&± are the boundary IR parameters. 
Recalling the bulk "lattice ^ IR" relation (see e.g. [53]) for the SSG model, 

A = ^ (4.14) 

and making the following identiflcation for the boundary "lattice -H- IR" relation for the 
BSSG model, 

V± = T^(s±(l + 2A)-2a±A) , (4.15) 

t9± = Avrfei. (4.16) 

one flnds the following. 

The above relations between the boundary IR and lattice parameters are among the main 
results of this paper. Such a relation was found in [S3] for the BSSG models with Dirichlet 



^See equation (4.35) of this reference. 

12 



boundary conditions. Note that (14.161) resembles the corresponding expression found for the 
BSG model in [38]. However, for the SG model, the bulk "lattice <-)■ IR" relation is given by 
A = -j^ instead of fl4.14p . Furthermore, we can also conclude that, combining the results 
(14. 6p . (14. 8 p and (I4.17p . that the NLIE generates the following SSG boundary S matrices 

R{9h ; A , r/+ , ^+) ~ P+(0, , -r/+ , ^+) P{9h) 

R{eh-A,yi-.^-)-P+{Oh.Tl-,^^-)P{Oh). (4.18) 

which agrees with the boundary S matrix proposed by Bajnok et al. in [50j for the BSSG+ 
model. This is also our main result. As pointed out earlier, the BSSG boundary 5* matrix 
is a product of SG and RSOS boundary S matrices. 

5 Discussion 

In this paper, we analyze the IR limit of a set of NLIEs that describe the BSSG model for a 
state of one hole. These NLIEs were obtained previously in [54] by taking the continuum limit 
of the NLIEs of the inhomogeneous open spin-1 XXZ quantum spin chain where the lattice 
boundary parameters (that form general integrable boundary terms as opposed to diagonal 
integrable boundary terms treated in [53]) satisfying a pair of real constraints (I2.10p . Thus, 
these NLIEs should correspond to more general boundary conditions than that considered 
in [53] . We computed the boundary S matrix which is the product of RSOS and SG terms. 
The SG term is then used to propose the relation between the boundary lattice and IR 
parameters. 

One could also further investigate this model. One could investigate the excited states of 
the model. It should also be possible to carry out such analysis for open spin-1 XXZ chain 
with non-diagonal boundary terms, the solutions to which are given in [16]. In contrast to 
the solution used in this paper, the solutions given in [16] are not restricted by any type 
of constraints among the lattice boundary parameters. It would be interesting to work out 
the NLIEs for these cases as well and analyze their UV and IR limits. Also, it would be 
interesting to see if there is an integrable spin-1 model that corresponds to the BSSG" model, 
explored in [5^. We hope to be able to address some of these issues in the future. 
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